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Rotational Brownian motion of colloidal magnetic particles in ferrofiuids under the influence of 
an oscillating external magnetic field is investigated. It is shown that for a suitable time dependence 
of the magnetic field, a noise induced rotation of the ferromagnetic particles due to rectification 
of thermal fluctuations takes place. Via viscous coupling, the associated angular momentum is 
transferred from the magnetic nano-particles to the carrier liquid and can then be measured as 
macroscopic torque on the fluid sample. A thorough theoretical analysis of the effect in terms of 
symmetry considerations, analytical approximations, and numerical solutions is given which is in 
accordance with recent experimental findings. 
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C3 . I. INTRODUCTION 
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CZ3 . 

Rotational Brownian motion of colloidal particles is a classical subject in statistical physics pj. Contrary to its 
translational cousin, i. e. "usual" Brownian motion, it does not result in changes of the particle location and is hence 
not as easily demonstrated in experiments. If, however, the suspended particles carry an electric or magnetic moment 
' their orientation couples to external fields and the electric or magnetic relaxation properties of the suspension are 
direct consequences of the rotational diffusion of the constituting particles. In particular, in the well developed theory 
of ferrofiuids the importance of rotational Brownian motion for the theoretical description of the often complex and 
O , surprising hydrodynamic and magnetic properties was recognized already a long time ago 0. 

Ferrofiuids are suspensions of nano-size ferromagnetic grains in carrier liquids like water or oil. They combine 
y—{ \ the hydrodynamic properties of Newtonian fluids with the magnetic behaviour of superparamagnets Q . Due to the 
^ ■ viscosity of the carrier liquid there is a coupling between the rotation of the magnetic grains and the local vorticity 
of the hydrodynamic flow. This coupling can, e.g., be used to spin up ferromagnetic drops floating in a non-magnetic 
On : fluid of the same density with the help of a rotating magnetic field 0, |{| . Complementary, a rotational ferrofluid flow 
exposed to a static magnetic field exhibits an enhanced shear viscosity @ . Various related effects such as "negative" 
rotational viscosity Q , magneto- vortical resonance |3| , and anomalously enhanced ac-response due to coherent particle 
rotation [t| have been investigated. They rely on the exchange of angular momentum between rotating particles and 
an oscillating magnetic field. In these cases, the imposed non-zero flow vorticity of the ferrofluid is crucial for breaking 
"jS ■ the symmetry between clockwise and counter-clockwise particle rotation. 

In the present paper we investigate a much more indirect and subtle aspect of the interplay between rotational 
Brownian motion of ferrofluid particles and their relaxation dynamics in an external magnetic field. We will show 
that a suitably designed time dependent external magnetic field without net rotating component may rectify the 
fluctuations of the particle orientation and set up a noise-induced rotation of the ferromagnetic grains. We will hence 
investigate how angular momentum can be transferred from an oscillating magnetic field to a ferrofluid at rest. The 
q ] effect was predicted theoretically and demonstrated experimentally in a previous short communication [ljj. I n t nc 
• • ■ present paper we provide much more details on the theoretical description and add several new results. 
. £h ' The extraction of directed motion from random fluctuations is an old and controversial problem in statistical 
^ . mechanics with a long and interesting history 0,0,0,0- Although excluded by the second law of thermodynamics 
for equilibrium systems, rectification of fluctuations is possible in systems driven sufficiently far away from thermal 
equilibrium |l5l Il6l | . The problem has gained renewed attention under the trademarks of "thermal ratchets" and 
"Brownian motors" due to its possible relevance for biological transport 0, and the prospects of nano-technology 

[UlUiil. 

Ferrofiuids are ide al sy stems to investigate such fluctuation driven transport phenomena and also to demonstrate 
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them experimentally |l(]|: As already discussed above, the rotational dynamics of the ferromagnetic grains is strongly 
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influenced by thermal fluctuations. Appropriate time-dependent potentials can be easily designed with the help 
of external magnetic fields. Finally, directed rotational transport in fcrrofluids should manifest itself as systematic 
rotation of the ferromagnetic nano-particles. This in turn can be easily detected from the resulting macroscopic torque 
on the carrier liquid. Various somewhat related but still quite different phenomena in rotational dynamical systems 
are treated in pi I2II2I ElE^ . 

Throughout the paper we will use two basic approximations which simplify the analysis considerably and which are 
rather common for ferrofluids The first is to neglect Neel- relaxation of the magnetization, i.e. the rotation of 

the magnetization vector with respect to the ferromagnetic particle. This is justified for particle sizes that are not too 
small and amounts to assuming that the magnetic moments are firmly attached to the geometry of the particles. Any 
reorientation of the magnetic moment hence requires a rotation of the particle as a whole. The second approximation 
is to neglect dipole-dipole interactions between the particles. Although these interactions may be important in 
concentrated ferrofluids, they are negligible in sufficiently diluted ones. The ratchet mechanism of central interest 
in the present investigation can operate without any particle-particle interaction. We will hence assume that our 
ferrofluids are sufficiently diluted and use a single particle model. For a more quantitative assessment of the role of 
the dipole-dipole interactions see e.g. |27^ and references therein. 

The paper is organized as follows. In section^ we introduce the general framework for describing the rotational 
Brownian motion of a single suspended ferromagnetic particle in an external magnetic field. Section lTTTl gives a detailed 
account of the symmetries characteristic of our system and discusses how these symmetries influence the possibility for 
a ratchet mechanism to operate. In section llVl the central Fokker-Planck equation describing our system is analyzed. 
Besides an investigation of the weak noise limit we describe the numerical solution of the Fokker-Planck equation 
and compare the results with two different approximate treatments yielding analytical expressions for the transferred 
angular momentum. Finally, section Ivl contains some conclusions. 



II. GENERAL FRAMEWORK 



We consider the rotational motion of a spherical particle of volume V and magnetic moment fi immersed in a liquid 
with dynamic viscosity rj and subject to a horizontal, time dependent, spatially homogeneous magnetic field H. The 
field is composed of a constant part H x parallel to the x-axis and an oscillatory part H y (t) with period 2tt/lu along 
the y-direction 

H = (H x ,H y (t),0) , H y {t + 2TT/u>) = H y {t) . (1) 

Different choices for the time dependence of H y are of interest. In the present paper we will mainly discuss two cases 
which are each representative for a whole class. Our first standard choice is 

H y (t) = H y 1} cos(wi) + H { y 2) sin(2wt + 5) , (2) 

(1 2) 

where the amplitudes H y ' , and the phase S are control parameters. The main features of this time dependence are 
a zero average over one period and the presence of a higher harmonic of the basic frequency. As a second example, 
we will also discuss the form 

H y {t) = Hf)+H^cos{ut), (3) 

for which the average over one period is different from zero. It arises naturally if the constant field component is not 
perpendicular to the time dependent one. In any case, the magnetic field is of a pure oscillatory character, i.e. it does 
not contain a net rotating component. 

The orientation of the particle at time t is described by the unit vector e(i) = /j,(t)/fi where fj, denotes the modulus 
of the magnetic moment. The time evolution of e is given by 

d t e = n x e , (4) 

where S~i(i) denotes the instantaneous angular velocity of the particle. Changes of fl are due to torques on the particle. 
Denoting by V the angular part of the three dimensional Nabla operator, the magnetic torque 

N mag = — e x VU = /i e x H (5) 

derives from the potential energy 



U(e,t) = -[i e • H(t) 



(6) 
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of a magnetic dipole in an external field [28j. Further, the viscosity rj of the carrier liquid gives rise to a viscous torque 

M 

N visc = -6r?Uft . (7) 

Additionally, the interaction between the rotating particle and the surrounding liquid also causes thermal fluctuations 
which generate a stochastic torque |2^ 



N stoch = y/UrjVkBT £(t) . (8) 

Here, is a vector of independent, (5-correlated Gaussian noise sources of zero mean, the noise intensity is related to 
the temperature T and the dissipation (J7J of the carrier liquid by the fluctuation-dissipation relation, and k B stands 
for Boltzmann's constant. Denoting the moment of inertia of the particle by /, the equation of motion for Jl acquires 
the form 



Id t fl + 6r)Vfl = n e x H + y/l2rjVk B T £{t) . (9) 

Eqs.QJ and form a closed set of equations for the description of the rotational motion of the particle. Using 
experimentally relevant parameter values (density of the particle p ~ 4 • 10 3 kg/m 3 , particle radius R ~ lOnm, 
viscosity -q ~ 1CU 1 . . . 10~ 3 Pas, time scale Sl/dt^i — 1CU 4 s) we find that the first term on the l.h.s. of eq. is five 
to seven orders of magnitude smaller than the second one. We may hence safely neglect inertial effects and find 
for fl in the overdampcd limit 

n = J^exH + V2DZ(t), (10) 
brjV 

where we have introduced D — k B T /6r)V . Using this result in Q yields a closed equation for the time evolution of e 

d t e = — (e x H) x e + V2D $(t) x e . (11) 
6r]V 

It is convenient to introduce dimensionless units. To this end we measure time in units of the inverse of the external 
driving frequency, t t/w, use QrjVuj/n as unit for the magnetic field strength, H > (6r]Vu)/fi) H, and rescale the 
noise intensity according to D i— > u> D. The evolution equation for the orientation e then reads 

d t e = (e x H) x e + V2D £(t) x e . (12) 
Introducing the Brownian relaxation time 

TB = Wf ' (13) 

we note that the rcscalcd, dimensionless noise intensity D occurring in (|12|l just gives the ratio between the relevant 
deterministic and stochastic time scales in the system: 

k B T 1 , , 

D = — — = . (14) 

GijVlu 2tb u> 

To proceed, we parametrize the orientation e of the magnetic particle by two angles 9 and ip according to 

e= (sinflcos^jSin^sin^jCOsfl) . (15) 
From l|12|l we then find the following Langevin equations for the time evolution of these angles [30| 

dt6 = -d e U + D cot 6 + V2D£ e (t) (16) 

1 y/2D 



d ^ = -^ U+ ^9^ t) - (17) 
In the dimensionless units adopted, the noise intensity D is given by l|14|) and the potential JBJ takes the form 

U(0,<p,t) = -sin (H x cos<p + H y {t) sin <p) . (18) 
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The thermal fluctuations £e(t) and £, v (t) are given by two independent, ^-correlated Gaussian noise sources of zero 
mean. Note that in (JT^ we are dealing with multiplicative white noise but since the multiplicative function is 
independent of (p(t), no ambiguity (Ito- vs. Stratonovich-interpretation j^]]]) arises. 

It is instructive also to consider the simplified situation of a one-dimensional dynamics in which the vector e is 
assumed to be constrained to the horizontal x-y-plane. Accordingly, l|16fl is replaced by 9 = ir/2 and (|17|l simplifies 
to 

d t <p = F{<p,t) + V2Dt v (t) (19) 

where 

F(ip,t) = -H x sin<^ + H y (t) cosip . (20) 

The observable of foremost interest in the present investigation is the time and ensemble averaged torque J^J exerted 
by the magnetic field upon the magnetic particle in the long time limit, i.e. after initial transient have died out. Since 
the magnetic field is constrained to the x-y plane, only the z-component of this magnetic torque can be different 
from zero. Suppressing the subscript mag in (J5J from now on, we denote the averaged z-component of the magnetic 
torque by (N z ), where (...) stands for the ensemble average over the different realizations of the noise terms in JT5J, 
(|17|l and the overbar represents the time average over one period of the magnetic field. Using Q and JSJ) we get 

1 



(N z ) = — I dt{shx9(t){-H x smtp(t) + H y {t) can(p(t))) . (21) 



o 



Exploiting l|18|) and i|17fl one readily finds the equivalent expressions 



2tt 2tt 

!/•.,„ 1 



(N z ) = ^Jdt (d<pU(6(t), <p(t),t)) =7^J dt ( Q Mt) si n 2 *(«)) (22) 
o o 

For reasons of ergodicity, the ensemble average in (|21|l is equivalent to a time average of a single realization over an 
infinite time interval. Then the extra time average over one period of the external driving drops out and we are left 
with 



t 

1 



(N z ) = lim — - — / dtdMt) sin 2 9(t) (23) 

(t f -ti)-*oo tf - ti 



ti 



and similarly for the equivalent expressions in Q2ip. (|22|l . 

In the absence of thermal fluctuations, the particle orientation is governed by the overdamped deterministic relax- 
ation dynamics given by I|17fl with D = 0. Hence, at any given moment, the orientation e(t) tends to align with 
the instantaneous magnetic field H(t) which lies in the x-y-plane. One can then show |32| that 9(t) converges to 7r/2 
for t — > oo and that (p(t) approaches a periodic long time behaviour. This implies via (|23|l that 

tp(t f ) - ip(U) 



(**)=, I™ V r v ' ; =0 (24) 

i.e., in the absence of thermal fluctuations no particle rotation will occur and no average torque can arise. An explicit 
example is displayed in Fig. ^ 

This scenario changes fundamentally in the presence of fluctuations, i.e. if D ^ 0. As shown qualitatively in Fig.^ 
for small noise intensities, the time dependence of <p(t) and 9(t) will still closely follow the deterministic trajectories, 
except that now also rare, fluctuation induced transitions between different deterministic solutions may occur. If the 
rates of forward (increasing tp by 2ir) and backward transitions (decreasing if by 2ir) are different from each other, 
then on the average a net rotation of the particle will occur, implying with 124fl that (N z ) ^ 0. This is a manifestation 
of the so-called ratchet effect [l5j in which an unbiased potential and undirected Brownian fluctuations cooperate to 
produce directed transport. The detailed operation of the ratchet mechanism for the rotational motion of a colloidal 
particle is the main focus of the present paper. 
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FIG. 1: Space-time plot of the potential t/(# = Tx/2,ip,t) from ea. 1181 for the time dependence © with H x — 0.3, .ffy 1 ^ = 
Hy 2 ^ = 1, and 5 = 0. Dark and bright regions correspond to small and large values of U, respectively. In the long-time limit, 
the deterministic dynamics 11611 . I17H with D = approaches 9(t) = n/2 and a periodical f(t), oscillating back and forth 
as represented by either of the full black lines. In the presence of a small amount of noise, occasional transitions across the 
unstable deterministic orbits shown as full white lines become possible which are schematically indicated by the dashed lines. 
The spatial asymmetry and temporal anharmonicity of the potential conspire to yield slightly different rates for noise induced 
increments and decrements of ip, respectively. As a result a noise driven rotation of the particles arises. For a detailed discussion 
see section llV Al 



III. SYMMETRIES 



Symmetry considerations admit conclusions from the equations of motion without actually solving them. For the 
investigation of ratchet systems, studies of symmetry turned out to be valuable since directed transport generic ally 
occurs if it is not forbidden by symmetries, a statement referred to as Curie's principle |l5l l33l l34l l35l 1361 1371 138. 39, 
liol fill |43| . Here we specifically address the question which transformations of the magnetic field components H x 
and H y (t) in leave the average z-component of the torque invariant, (N z ) 1— > (N z ), and which lead to its inversion, 

From Eq. (|23|l one readily concludes that the following transformations of tp and 9 imply the indicated transfor- 
mation behavior of (N z ): 





-> tp(t) + Aip 


and 




-0(t) 


implies 






¥>(*)•- 


-> -ip(t) + Aip 


and 


9(t), 


+ 0(t) 


implies 




- 






and 




-> -9(t)+n 


implies 








-> ip(t + At) 


and 


9(t)H 


-> 9(t + At) 


implies 




- (N z ) 




- ¥>(-*) 


and 


9(t), 


+ 0(-t) 


implies 







(25) 
(26) 
(27) 
(28) 
(29) 

where Ay and At are arbitrary. Observing that the values of both the original and the transformed 9(t) are restricted 
to the interval [0, it] and that the entire interval will indeed be sampled by some realization of the dynamics l|16|) . i|17[l . 
it follows that no transformations of the form d(t) i-> 9(t) + A9 with AO ^ or 9(t) i-> -0(t) + A6> with A9 ^ it are 
possible. Finally, in order to conclude (N z ) t— > — {N z ) in l|2"§|) one has to exploit that in (|2~3l) the limit (£/ — ij) — > 00 
is not necessarily tantamount to t/ — ► 00 with kept fixed, but can also be realized by letting tj — ► — 00 with kept 
fixed. 

Closer inspection of (|23|l reveals that any further symmetry transformation that leaves (N z ) invariant or changes 
its sign can be composed of the "elementary" transformations <|25|) - (|29ll . 

Next we ask ourselves what equations of motion govern the transformed 6(t), <p{t) defined by (|25|) - (|29|l . In the 
simplest case (|25|l one readily concludes from (|16|) and l|17(l that the transformed 9(t), (pit) satisfy the same equations 
of motion l|ltj|) and H17|) but with a transformed potential U(9, ip,t) 1— > U (9, ip — Aip, t) + AU with arbitrary A/7. 
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Turning to the second transformation (|26() one finds similarly that the transformed 9(t) and ip(t) satisfy the equations 
of motion (|16[) and l|17l) but now with a transformed potential U(9,ip,t) i— ► U(9,—p + A,t) + AU and in addition 
with a transformed noise i— > Since the statistical properties of the two noises — £^(i) and £ v (t) are equal, 

i.e. they are identical stochastic processes, the same follows for the corresponding processes given by H16[) and H17[) . 
In other words, to each realization of the original 9(t), cp(t) corresponds a realization of the transformed 0(t), cp(t) 
occurring with the same probability, and vice versa. Since the average torque in 1(23|1 is independent of the specific 
realization 9(t), ip{t) with probability 1 (i.e. up to a set of realizations of measure zero) we can conclude that in this 
case (N z ) i— > —(N z ). Similar considerations can be done for the cases given by (|27[) and l|28|) . The last transformation 
(|29[1 is insofar special as a mapping between the dynamics of the original and the transformed variables is possible 
only, if the term D cot 8(t) in l|16[) identically vanishes, i.e. if 6(t) = tt/2. We therefore find the following implications 
of transformations in the potential U(9,(p,t): 





-> U(9, ip - Acp,t) + AU 


implies 




->(Nz) 




(30) 


U(0,<p,t)* 


■+ U(6,-(p + A(f,t) + AU 


implies 








(31) 


U(0,(p,t)^ 


-> U (—9 + 7r, ip, t) + AU 


implies 








(32) 


U{B,<p,t)^ 


-> U(6,ip,t + At) + AU 


implies 








(33) 




+ -U(6,tp,-t) + AU 


implies 






(provided 6{t) = tt/2) , 


(34) 



where AU is an arbitrary constant. 

Finally, we have to determine all those transformations of the magnetic field components H x and H y (t) in JIJ which 
correspond via i|18|) to transformations of the potential U(9,<p,t) given in l|30l) - (|34[l . Since changes in H x and H y (t) 
always imply changes in the potential U that depend on 9 and ip, we can only induce transformations with AU = 0. 
The only non-trivial way to realize (|30|l is (H x , H y (t)) i— » (—H x , —H y (t)) combined with Aip = (2n + l)ir with 
an arbitrary integer n. Similarly, (|31|l and (|33|) are equivalent to (H x , H y (t)) i— > (H x , —H y (t)) and (H x , H y (t)) i— > 
(if x , + At)). On the other hand 13211 only admits the trivial realization (H x , H y (t)) i— > (ii^, H y (t)). Finally, 
l(34|l can be implemented by (H x , H y (t)) i— > (—H x , —H v (—t)). In summary, we find that 



(£T X , ff v (t)) h 


(-if x , -#„(*)) 


implies 




- w 




(35) 


(# x , «■„(*)) H 


+ (if,, -fl„(t)) 


implies 




- - w 




(36) 


(# x , ff y (t)) k 


- (ff x , H y (t)) 


implies 




- w 




(37) 


(# x , ff„(t))'- 


■+ (ff x , H y (t + At)) 


implies 




- w 




(38) 


(# x , H y (t)) H 


-> (-ff x> -#„(-«)) 


implies 




- -<JV.) 


(provided 9(t) = tt/2) , 


(39) 



where the trivial result l|37J) is listed only for the sake of completeness. 

While H35 p -Q38 p are intuitively more or less obvious, l|39|) respectively (|34|) is not. The additional condition 9{t) = tt/2 
shows that the constrained dynamics (|19|l . (|20|> has an extra symmetry which is lost if the full two-dimensional 
dynamics i|16|) . I|17|l is considered. 

We note that H35[) - l|39|) could also have been obtained in a more direct way. Our somewhat more involved line 
of reasoning has the advantage that we can exclude that there are any transformations other than 1)35(1 - (|39|l and 
combinations thereof which would leave (N z ) invariant or change its sign. 

In order to see that the above considerations already allow some non-trivial conclusions about the possibility of 
angular momentum transfer from the oscillating magnetic field to the particle we combine (|35l) and l|36[) . (|36|l and 
and and ((23) to obtain 



(H x , 


Hy{t)) - 


* {-H x , H y (t)) 


implies 




- -(N z ) 




(40) 


(H x , 


Hy{t)) H 


■+ (H x , -H v {t + At)) 


implies 








(41) 


(H x , 


H y {t)) H 




implies 




- - W 


(provided 9{t) = tt/2) 


(42) 



respectively. From (|40|) we immediately can infer that the constant field in ^-direction is indispensable for a non-zero 
average torque, 

H x = => jN7j = 0. (43) 

Similarly (|41|l implies that for a time dependence obeying H y (t) = —H y (t + At) for some At the average torque has 
to vanish as well 

H y (t) = -H y {t + At) JN7)=0. (44) 
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Choosing in particular At — tt, i.e. half the period of the external driving, we realize that an oscillating field H y (t) 
with a Fourier expansion containing only odd multiples of the basic frequency 

H y (t)= J2 cos(nt + S n ) (45) 

n=l,3,5,... 

will result in a zero average torque (N z ), irrespective of the particular choices of the amplitudes H y n ^ and phases S n . 
Eas-H43p and l|45|l motivate our choices J2J and J2J as simple time dependencies of the oscillating field resulting in a 
non-zero average torque (N x ). 

For the special time dependence © we find from (|41ll 

JN^(S + tt) = -Jn7)(S) (46) 

As a consequence, upon continuously varying S we can infer that there must exist a, Sq € [0, tt) such that (N z ) = 
for S = Sq or 6 = So + tt. It is important to note that this case of zero averaged torque is qualitatively different from 
the situations described by eqs. I|43|l and (|44|l since it is due to the fine tuning of a parameter rather than resulting 
from an underlying symmetry. Consequently, upon variation of S around So the average torque (N z ) changes sign, 
which is a realization of a so-called current inversion 15]. Furthermore, by fixing S = Sq, a sign change of (N z ) will 
generically also occur upon variation of any other parameter of the system. Interesting examples for such a parameter 
include the driving frequency u>, the noise intensity D, and the particle volume V. In the latter case we hence face the 
intriguing possibility that in a polydisperse ferrofluid under the same experimental conditions the larger and smaller 
particles will rotate in opposite directions. 

We finally note that in the case where the dynamics is constrained to the x-y-axis, the additional symmetry 
(|4*2|) implies that So = n/2. Consequently the S- values where the torque changes sign are fixed by the additional 
symmetry and do not depend on the other parameters of the problem. Hence no current reversal upon changes in 
these parameters is possible. This is the main difference between the full two-dimensional dynamics I|17|) and 

its simplified one-dimensional counterpart 1|19[) . see also [37tl44l l45l |46| . 



IV. ANALYSIS OF THE FOKKER-PLANCK EQUATION 

A quantitative analysis of the overdamped Brownian rotation of a single ferrofluid particle in a time dependent 
external field can be built on the numerical simulation of the Langevin equations (|16l) . 1|17[) . Such simulations 
unambiguously show the occurrence of noise induced rotation of the particle in an oscillating field [3(j . For a detailed 
study of the influence of the various parameters of the problem on the transferred torque (N z ) it is, however, more 
convenient to start with the equivalent Fokker-Planck Equation (FPE) for the probability density P = P(9, ip, t) for 
the orientation l|15|l of the particle 

d t P =V(PW) + L>V 2 P 

= 4— d e (sin 9Pd g U) + d v (P d v U) + — a dg (sin 9d g P) + 3%P . (47) 

sin sin sm9 S m v 

Here V denotes the angular part of the three dimensional Nabla operator. The equivalence of Ijlfijl . (|17|l and (|47|l is, 
e.g., shown in .30|. For a time-periodic potential U(9,ip,t), the solution P(0,ip,t) of (|47|l will also be periodic after 
initial transients have died out. From this solution we can determine the average orientation of the particle 

/*27T rn 

(e(t)) = / dip d9 sin #(sin6> cos tp, sin 9 sin ip, cos 0) P(0,ip,t) (48) 
Jo Jo 

with the help of which the average torque from (|21|I - 1|22J) can be calculated according to (0) in yet another way, namely 

W*) = l^J*dt[{e(t))xlI(t)} z . (49) 

In the present section we discuss the numerical and approximate analytical solution of the Fokker-Planck equation 
(|4TJl and the corresponding results (|4*9"|) for the averaged torque (N z ). To first substantiate the intuitive un- 

derstanding of the ratchet effect in our system advocated with Fig. ^ we start by considering the weak noise limit 
D — > 0. To keep the analysis simple we will restrict ourselves in this part to the simplified one-dimensional model 
defined by i|19|) , (|2(J|) . In the second subsection we detail the numerical methods used to solve the FPE for the general 
case. Finally we present some analytical results from the perturbative solution of the FPE. 
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A. Instantons in d = 1 

Our starting point is the simplified one-dimensional model dynamics <|19|) . 1|20|) with a ^-correlated, unbiased Gaus- 
sian noise source £, v {t). The deterministic dynamics described by (|19fl with D = has a family of stable and unstable 
periodic orbits f s (t) ± 2nn and <p ns (t) ± 27rn, respectively, where n is an arbitrary integer (cf. Fig. ^) 32|. Hence 
deterministically the particle simply follows the direction of the magnetic field with a certain time lag due to viscous 
damping and no particle rotation can occur. 

For any non-zero noise intensity D > 0, there will be few stochastic transitions between the deterministic trajectories 
and their detailed form in the limit D — > can be determined. Comparing the associated rates for forward (increasing 
ip by 2tt) and backward (decreasing ip by 2ir) phase slips we may directly obtain direction and magnitude of the 
noise-driven transport by means of l|24|l. 

The analysis of this section builds on the path integral representation for the transition probability of a Markovian 
stochastic process |47l l48l 149^ . The main complication is the time dependence of the potential U(8,<p,t). In notation 
and general strategy we follow the recent detailed analysis of noise-driven escape over oscillating barriers [50j . We will 
only determine the dominating exponential term in the transition probabilities. A more detailed analysis including 
also the prefactor can be done along the lines of ref. [5(| • 

The transition probability p(<Pf,tf\cpi,ti) of the Markov process described by (|19|l can be written in the form 

<p(*f)=<Pt 

p(Vf,t f \cp i ,t i )= J ^(-)e X p(-^P) , (50) 

¥>(t»)=Vi 

where the action functional S" [<£>(•)] is given by 

%(•)] = \ £ dt(dMt) - F(<p(t), t))\ (51) 

with F(ip, t) defined in H20JI . The intuitive understanding behind this representation is that all possible trajectories (p(t) 
starting at initial time U in ipi and arriving at final time tf in ipf contribute to the transition probability p(<Pf, tf\(pi,ti), 
each with a weight related to the value of the action S evaluated along this trajectory. In the weak noise limit D — > 
a Laplace argument tells us that only those transitions <p*(t) contribute significantly to p(ipf,tf\tpi,ti) for which the 
action is minimal. All others have probabilities which are exponentially small for D — ► 0. In the weak noise limit 
we are hence able to determine the precise form of the (dominating) stochastic transitions by solving a variational 
problem. 

Under rather general conditions which are satisfied in our particular setting one can show that the transition rate 
between two stable orbits of the dynamics is determined by the probability to reach from the starting stable orbit the 
unstable orbit separating the two stable ones |5(j ■ The subsequent relaxation to the final stable orbit is then purely 
deterministic and does not contribute to the transition probability p (<p>f-, tf\<pi, U) in the weak-noise limit. 

We have hence to solve a variational problem for which the initial and final point of the trajectory are not fixed. 
Instead they are required to lie on two known functions (f s (t) and tp ns (t), respectively. If we perform the optimization 
with respect to these locations we find 

d t <p*(t i ) = F(<p*(t i ),t i ) and d t sp*{tf) — F(ip*(tf),tf). (52) 

Hence the optimal transition trajectory <p*(t) must be parallel to the stable and unstable orbit at the contact points. 
This in turn implies U = — oo and tf = oo. The optimal transition trajectory if*(t) hence starts out at i, = — oo at 
the stable orbit, moves for a very long time near to it, changes then in a rather short time interval to the immediate 
neighbourhood of the unstable orbit (hence the name "instanton" ) which it finally reaches at tf — oo 50] . 

From the minimization of the action #[<£(•)] with respect to the functions (p(t) we get the usual Euler-Lagrange 
equation for the minimizing instanton (p* (t) which for the action given by l|51|) takes the form 

%<p*(t)=a t F(<p*(t),t)+F(v*(t) t t)d v F( V >*(t),t) . (53) 

Choosing for the time dependence of H y (t) the example specified by @ we have solved (|53|) numerically as a 
boundary value problem for a system of ordinary differential equations. The initial and final points were chosen 
slightly off the stable and unstable orbit, respectively. If the difference between the initial and final time, tf — t», is 
larger than several periods of the external driving, the solution hardly depends on the precise value of these deviations 
from the deterministic orbits. Having obtained tp*(t), the corresponding value of the action S can be calculated from 

(EU- 
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FIG. 2: Numericaly determined instanton trajectories (full lines) describing the dominating weak noise stochastic transitions 
between stable (dashed lines) and unstable (dashed-dotted lines) deterministic periodic orbits of the one-dimensional rotational 
dynamics of a ferrofluid particle described by 119H . 120H . The parameter values are H x = 0.3, H^ l) = = 1 and S = 0. 
For the forward instanton, increasing ip by 2n (right), the value of the action is S = 0.21472, for the backward instanton, 
decreasing ip by 2n (left), it is S = 0.30969. Hence, a noise induced average net rotation in the "forward direction" results in 
the weak noise limit, implying (N z ) > according to 1241 1. 



In fig- Elthe resulting instanton trajectories p*(t) corresponding to backward and forward transitions are displayed. 
They show the qualitative behaviour discussed above. Note, however, that only the noise induced transition from 
the stable to the unstable orbit is shown. The complete transitions giving rise to ip 1— > ip ± 2ir comprises also the 
subsequent relaxations from the unstable orbit to the next stable one. Consequently the stochastic transitions show 
substantially more structure than the rather poor caricatures used to represent them in Fig. ^ 

The corresponding values of the action are S — 0.21472 for the forward instanton increasing ip by 27r and S — 0.30969 
for the backward instanton decreasing <p by 27r. Hence, at least in the weak noise limit, forward transitions are more 
likely (their action is smaller) than backward transitions and correspondingly for the parameter set specified in the 
caption of Fig. El we find a noise driven rotation in the direction of increasing ip. The mismatch between forward and 
backward transition clearly demonstrates the operation of the ratchet effect in our system. With the relevant transition 
probabilities being of order exp(— 0(1) /D) the effect is however very weak and in particular rapidly disappears with 
D^0. 



B. Numerical solution of the FPE 



For the numerical solution of the FPE describing the complete two-dimensional dynamics for general values of the 
parameters it is convenient to expand the probability distribution P(0,ip,t) in spherical harmonics Y 1 m (9,ip) |5lj 

oc I 

P(8, = £ £ A ^{t) Yr(6, <p) , (54) 

2=0 m=— i 

with so far unknown time dependent complex expansion coefficients A^ m (t). Since P(9,ip,t) is real we have 

A- m (t) = (-l) m Al m (t) , (55) 

with the star denoting complex conjugation. Moreover, from the normalization of P(0, <p>, t) we find 

Ao, (t) - -J= . (56) 

V47T 
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Plugging the ansatz l|54|) into the FPE (|47(l , exploiting (|18fl , and using well-known properties of the spherical harmonics 
we find for I > 1 the following set of ordinary differential equations 

d t A lim (t) = - Hx (h^l - l, m - 1) An, m -!(() + h 2 (l + 1, m + 1) A i+1 , m _!(*)) 



ff a + zff„(t) 
2 

Dl(l + l)A„ n (t) (57) 



(hi (I - 1,-m - 1) A/_i, m+ i(f) + /i 2 (Z + 1,-m- 1) A m , m+ i(t)) 



where we have introduced the auxiliary factors 



h n \ fij.<>\ (l + m + l)(l + m + 2) 

hl (l,m) = (/ + 2) J (2/ + 1)(2/ + 3) (58) 



, \ (I — m)(l — m — 1) 



Moreover, we find for the ensemble averages of the components of e(i) from l|15|l 

(e x (i)) - (sin0(t) cos^(i)) = -J*jL%A ltl (t) (60) 
(e,(t)) = (sin0(*) sin^t)) = ^9fAi,i(t) (61) 

(e,(t)) = (cos 0(i)) = ^ A ll0 (t) (62) 

where 5iz and Q'z denote the real and imaginary part of the complex number z, respectively. For the z-component of 
the ensemble averaged but still time-dependent magnetic torque (N z ) this implies via (|5J) 




H y {t)nA hl {t) + H x $A hl (t)) . (63) 



We have used two somewhat different procedures to explicitly calculate this torque. The first is more general, the 
second is slightly more efficient for the special case of the time dependence ^ . 

In the first method we solve the system (|57H of ordinary differential equations for the expansion coefficients A; )Tn (i) 
with I < Z max by a standard routine. We have used values of Z max up to 25 but found typically 5 < / max < 10 to be 
sufficient. Starting with some initial condition, we integrate the equation until the periodic time dependence of the 
solution is reached, which was usually the case after at most two periods of the external field. We then integrate further 
for one period of the external driving, use the result for Ai t i(t) in (|63l) . and calculate the remaining time average 
numerically. Some care is needed here since this time- averaged torque (-/V z ) is typically 2-3 orders of magnitude smaller 
than the typical time-dependent values of N(t). In the time average we are hence subtracting numbers of equal size 
from each other giving rise to well-known problems of numerical accuracy. The advantage of this method is that it 
works for any numerically sensible time dependence of the magnetic field H y (t) in Q • 

The second method builds on the fact that the long-time solution of the FPE will be periodic in time with the 
period 27r of the external driving. It is hence useful to expand the time-dependent expansion coefficients Ai^ m (t) in 
(|54"Jl into Fourier modes with respect to time 

oo 

Ai, m (t)= A i,m, s e ist - (64) 

s— — oo 

Instead of a system of ODE's for the coefficients A^ m (t) we find now from the FPE l|47|l a system of algebraic equations 
for the coefficients Ai irlltS . Moreover, the time averaged torque (N z ) can be readily expressed as function of some 
of these coefficients by means of (JS3J. The advantages of this method are that effectively we are directly dealing 
with the stationary solution of the FPE, i.e. no initial conditions are necessary and no equilibration process must be 
simulated, and that the time average of the torque needs not to be performed numerically. The disadvantage is that 
we have to specify H y (t) in order to find the algebraic system for the Aj m s and that for general H y (t) this system 
will be rather dense with respect to the index s. 
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FIG. 3: Left: Averaged torque (N z ) for the time modulation @ with H x — 0.3, Hy 1 * 1 — H y 2 ^ = 1 and D — 0.2 as function of 
the phase angle 5 as obtained from the numerical solution of the FPE (full line). Symbols are results from simulations of the 
Langevin equations 1161 and (11711 as reported in [30j with statistical errors being smaller than the symbol size. The squares 
on the horizontal axis are the points (n/2,0) and (3n/2,0) at which the current reversal would occur in the reduced model 
IT51 . Right: Depend ence of (N z ) on the relative particle size for the same values of the parameters and S fixed to r5o — 1.73, 
the first point of current inversion in the left panel. The full line corresponds to a proportional change of the magnetic and 
hydrodynamic radii (top sketch), the dashed line is for a change of the hydrodynamic radius only (bottom sketch). 

For the special time dependence the situation is somewhat more gratifying since H y (t) involves only two Fourier 
components. Hence there are only couplings between coefficients A^ m ^ s differing in their s-value by at most 2 and the 
result for (N z ) is a linear combination of just 5 different coefficients Ai tTn ^ s . This makes the numerical analysis rather 
fast. Both methods must of course yield the same results when applied to the same H y (t). We have frequently used 
both to verify our numerical findings. 

Using the numerical methods described above we have verified all the predictions derived in secti on ITTll on the 
basis of symmetry arguments. In particular we find (N z ) = if H x = in accordance with (|43fl and (N x ) = for 

H y (t) = cos(t), i.e. for 10) with H y 2) = 0, or H y (t) = H y 1} cos(i) + H ( y 2) cos(3i) in accordance with (|4*Ijl . 

In Fig. [3] we give some examples for current inversions in our system. The left panel shows the averaged torque 
(N z ) as function of the phase angle S for a time dependent magnetic field H y {t) of the form (J2J. The results are 
consistent with Moreover, the values of 5 at which the torque changes sign are clearly different from ir/2 and 

3ir/2 indicated by the squares on the horizontal axis. Hence the numerical solution verifies the statement of section UTTI 
that the current reversal in the full two dimensional model i|ltj|) , (|17f) occurs at values of S which are not fixed by any 
symmetry but rather depend in a complicated manner on all the other parameters of the problem. 

This is further corroborated by the right panel of Fig. [31 depicting the torque (N z ) as function of the relative 
particle size for S = do- There are two different ways to scale the particle size that make sense in a ferrofluid. In the 
first, both the ferromagnetic core of the particle as well as the polymer coating are scaled with the same factor. This 
gives rise to a proportional change of the magnetic moment /i and the particle volume V. In the dimensionless units 
adopted this implies unchanged fields H x and H y (t) and a noise intensity D rescaled according to l|14[) . In the other 
case, only the polymer coating is scaled and therefore only the hydrodynamic radius of the particle changes. Hence 
the magnetic moment remains the same whereas the volume V changes. Correspondingly the dimensionless fields H x 
and H y (t) and the noise intensity D change by the same factor. In view of (|47|l this can be absorbed in a rescaled 
time which in turn is equivalent to a changed frequency u of the external field. As anticipated in section ITTT1 we find 
current inversions in both cases. 



C. Effective field approximation 



It is very helpful to have some analytical expression for the average torque (N z ) since the dependencies on the various 
parameters of the system can then be identified much more directly. In the present and subsequent subsections we 
discuss two approximate methods to calculate the time averaged torque in our system. 

The first one employs the so-called effective field method to approximately solve the FPE. It has been a standard 
tool in the theory of ferrofluids for many years |52(. Two main ingredients are necessary. First, let us recall that the 
stationary (equilibrium) solution of the FPE l|47(l for a time independent homogeneous magnetic field H is of the form 
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where H denotes the modulus of the magnetic field, H = |H|. Averages with this distribution are easily calculated, 
in particular we find the well-known result 

(e)(°) = L(H/D) H/H (66) 

with the Langevin function 

L{x) = coth(a;) - 1/x . (67) 

Second, for general, time-dependent fields H(i) the FPE ljT7|) can be used to derive the exact equation 

d t (e) + 2D(e) = -(e x (e x H)) (68) 

for the time evolution of the average (e). Note that this equation for the first moment of P(e,t) has the usual flaw 
of being not closed but involving higher moments of the probability distribution. 

The central idea of the effective field approximation is now to assume that the solution P(e, t; H(i)) of the FPE for 
a general time dependent magnetic field H(i) can be written as equilibrium distribution for some a-priori unknown 
effective field H e (t), i.e. P(e,i;H(t)) = P' ' (e; H e (t)). The average on the r.h.s. of II68II can then be performed and 
this equation gives rise to an evolution equation for the effective field H e of the form |52| 

a , { mm h„) = _ 2B £ffiZ£i (Ht - h, - ".-y/"> (Ht , |H , x H » , 69) 

Given H(i), this is a closed equation for the time evolution of H e (£). Having obtained H e (t) we get the desired 
approximate result for (e) from Ij66|l 

(e) e = L(H e /D)H e /H e . (70) 

The solution of the partial differential equation l|47() is hence replaced by the solution of the set of three coupled 
ordinary differential equations l|69|) . There are other, equivalent forms of l|fc>9|l . see e.g. [S^. Exploiting l|69|) . it is also 
possible to derive the following equivalent closed equation for (e) 0] 

d t (e) e = -(1 - ^)(H e - H) + (1 - 3D^)(|^ • (H e - H)) @± , (71) 

where (e) e denotes the modulus of (e) e and the effective field H e is to be expressed as function of (e) e by the inverse 
of 17U|) . Note the difference between (e) e = |(e) e | and (|e|) e = 1. 

In general, (|69|l can only be solved numerically. In Fig. 0] we give a comparison between results obtained in this 
way for the time dependence @ and the corresponding outcome of the numerical solution of the FPE. It is clearly 
seen that the effective field approximation yields rather accurate results for the time dependent orientation of the 
particle and, equivalently, for the time dependent magnetization of the ferrofluid. On the other hand the results for 
the time averaged torque (N z ) differ from the numerically exact values by a factor between 2 and 3. The reason for 
this discrepancy lies in the fact that (N z ) is typically 2-3 orders of magnitude smaller than typical values of the time 
dependent torque (N z ). In calculating the time average of the torque we hence subtract quantities of comparable 
order and so amplify the inaccuracies of the effective field approximation from a few percent to a few hundred percent. 



In order to obtain an analytical expression for (N z ) within the framework of the effective field approximation, it is 
useful to consider the case of small values of H/D in which the Langevin function L(x) or equivalently its inverse may 
be approximated by the first terms of the respective power series expansion. Using the lowest order approximation 
L(x) ~ x/3 we have H e ~ 3D(e) e and from l|71(l it then follows that 



For large t we therefore find 



2H X 
3D 



9 i (e) e + 2P(e) e = 2H/3. (72) 

(73) 



2H ( y ] 2D cos t + sin t H { y ] D sin(2i + 5) - cos(2i + S) 

~~3 4D 2 + 1 D 2 + l ' ( ' 
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FIG. 4: Horizontal components of the average orientation (e) of the particle as function of time for one period of the external 
magnetic field. The lines show the results for the x-and y-component of (e) (full and dashed line, respectively) from the 
numerical solution of the Fokker-Planck equation 147H with a time dependence of the external field as given in 0. The symbols 
are the corresponding results from the effective field theory. The parameter values are H x = 0.3, H y x ^ — H y 2 ^ = 1, 5 = as 
well as D — 1 (left) and D = 0.2 (right). The corresponding values for the time averaged torque (N z ) are also displayed. 



Hence (e x ) e tends asymptotically to a constant whereas (e y ) e contains only oscillating parts. From © we may thus 
infer that the torque (N z ) harmonically oscillates and therefore its time average vanishes, (N z ) = 0. 

To get a non-trivial result for (N z ) we have to push the expansion of the Langevin function further. In other words, 
the non-linearity of the magnetization curve is essential to get a non-zero average torque. If the next term in the 
expansion is included, L(x) ~ x/3 — x 3 /45, we already obtain an approximation for L(x) that is accurate to within 
one per cent for the experimentally relevant values of x which are generally less than 1. Calculating the corrections 
to H73fl and (|74fl induced by these higher order terms we finally find after some algebra the following result for the 
average torque 

— - = H x (hPYH™ cos5 + 2DsmS 

[ z)e 30 (1 + L> 2 )(1+4L> 2 ) 2 ' 1 ' 

This result is in accordance with our symmetry considerations of section ITTT1 since it shows that both a non-zero field 
component in x-direction and an even higher harmonic in the time dependence of H y (t) are essential for a non-zero 
torque to occur, cf. i|43[l . (|44|l . 



D. Perturbative solution of the FPE 



In the last subsection we saw that the effective field approximation alone is not sufficient to get an explicit analytical 
expression for the time averaged torque (N z ) in our system. Additionally we had to expand the Langevin function in 
its argument H/D which amounts to a small field or equivalently large noise expansion. Moreover, the approximate 
result for the torque was typically at variance with the numerically exact result by a factor of about 2. 

It is therefore tempting to avoid the effective field approximation altogether and to use a perturbation expansion in 
the ratio of magnetic field strength and noise intensity right in the Fokker-Planck equation 1J47J1 • This is conveniently 
accomplished by making the formal substitution H i— ► eH in H47[) and by expanding the solution of the FPE in a 
power series in e 

P(9, tp, t) = (0, ip,t)+e P (1) (6, ip, t) + e 2 (6, if, t) + e 3 P (3 *> (6,<p,t) + ... . (76) 

At the end of the calculation e is put equal to one. It w ill tur n out that we have to use the expansion up to third 
order in e to get a non-zero result for the average torque (N z ). Since the potential U(9,tp,t) is of order e we find in 
zeroth order 

li m p(°)(^^i) = i , (77) 
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describing the stationary distribution of pure rotational diffusion. Using the expansion (|76|l in i|47|) we get a system 
of linear inhomogcncous partial differential equations for the various pM(0,<p,t) that can be solved in principle one 
after the other. 

This procedure is, however, rather cumbersome. It can be somewhat simplified by using the exact equation (JfJSJ), 
the r.h.s of which is of order e. Hence to find (e) from this equation to order e™ it is sufficient to calculate the average 
on the r.h.s. with an expression for P(9, ip, t) accurate to order e*™ -1 -'. In our case the lowest order non-zero result for 
(N z ) is 0(e 4 ), hence we need (e) to order e 3 (cf. JSJ), (|15fl ) and therefore the second order result for P(8,ip,t) would 
be sufficient. 

It is, however, more convenient to use again the expansion (|54|) of P(6, tp, t) in terms of spherical harmonics. The 
perturbation ansatz (|76|l then translates into an expansion of the coefficients Ai iTn (t) in powers of e 

A hm {t) = Afi(t) + eA^(t) + e 2 A™(t) + e 3 A®® + ... (78) 

The ODE's (|57|l couple only coefficients Ai <m (t) which differ in both / and m by 1. Moreover, we know from (|56|) 
Ao.o(t) = 1/VAtt = 0(1) and hence A^ m (t) = 0(e l ). As shown by 163(1 the torque is determined by A\_i(t) alone. We 
hence need an expression for A\,i(t) correct up to order e 3 . 

To this end we determine from l|57|l and l|56l) first the 0(e) terms of Ai :in (t). With these we calculate the necessary 
A2, m (t) to order e 2 which in turn are used to determine the 0(e 3 ) terms in Ai t \(t). The procedure is sketched in 
n g- El Note that it is sufficient to calculate coefficients A^ m (t) with m > because of the symmetry property (|55|l . 
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FIG. 5: Diagrammatic sketch of the perturbative determination of the expansion coefficients Ai, m in 1541 . As shown by the 
arrows only Ai )m differing in both / and m by 1 are coupled by equations 1571 . In a first step we start from Ao,o as given in 
156> and obtain the 0{e) terms of A\ t ±\ (dashed arrows). Next ^2,0 and A2,±2 are determined to order e 2 (dotted arrows). 
Finally these are used to get the 0(e 3 ) contributions to (full arrows) from which the torque can be extracted to order e 4 
using I63H . All other Ai >m are irrelevant for (N z ) at this order in e. Note also that the left half of the diagram is redundant 
due to the symmetry property 1551 . 

In this way we first get from the 0(e) equations 

(0 t + 2D)A { ^(t) = -J=( Hx - iH y (t)) 
(d t + 2D)A§(t) = 0. 

Since we are interested in the asymptotic solution valid for large t we can put A^l(t) = 0. Next we use the 0(e 2 ) 
equations 

(d t + 6D)A^(t) 

(d t + ei>)43(t) 



-—(H x -iH v (t))A$(t) 



= \l~[H x MAV(t)-H y (t)ZA?>(t) 



1W1 
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FIG. 6: Left: Time averaged torque {N z ) as a function of the noise intensity D for the time dependence of the magnetic 
field. The parameter values are H x = 0.5D, Hy = 0.5D, H y 2 ^ = 0.15D, and 5 = 0. The scaling of the magnetic field strength 
with D ensures that the perturbation parameter e stays constant. Shown are the perturbative results 1801 (circles) and 17511 
(squares) together with the numerically exact result (full line). Right: Same for the time dependence of the oscillating field 
with -ffy ' = Hy 1 ^ = 0.5D. The approximate results are given by 1821 and I81H . respectively. 



(3) 

to determine the relevant A2, m {t) to the desired order. Finally we obtain A\ {(t) from 

(d t + 2D)A<?l{t) = -^{H x - iH y (t))A^(t) + ^{H X + iH y (t))A^(t) . (79) 

The explicit expressions for the various coefficients A are rather long, the above system of equations is therefore 
conveniently solved with the help of a computer algebra. Focusing on the specific time dependent field H y (t) from , 
exploiting (|(j3fl and performing the time average over one period of the external driving we finally get the following 
lowest order perturbative result for the average torque 

T W\2 17(2) 



H x (Hy ) 2 H y ' 9(1 + 29L> 2 ) cos 6 + 2D(1 + 99 D 2 ) sin<5 



40 (1 + L> 2 )(1+4L> 2 )(1 + 9D 2 )(1 + 36L> 2 ) ' ^ 

The dependencies on the parameters H x , Hip , and H^ are identical to those of the effective field result (|75J) . whereas 
the dependence on the noise strength D is more complicated. One readily verifies that the symmetry properties lj35|l - 
(|38JI are satisfied by H8()(l . In particular, these symmetries imply that the torque (N z ) must to be an odd function of 
both the static field and the oscillating field. Moreover, symmetry reasons also imply fl5l | that in linear order of the 
oscillating field one cannot expect a finite (N z ), Rather, a coupling of several modes of the oscillating field (harmonic 
mixing) is required, see also [33 144|, |46j . All together, these arguments explain that the first nontrivial contribution 
must be (at least) of fourth order in the magnetic field strength and that within a linear response theory one will 
always find (N z ) = On the other hand this does of course not imply that the non-linear magnetization curve 
alone would "explain" the effect. To understand the microscopic origin of the angular momentum transfer from the 
external field to the particle one really has to go to the description in terms of the FPE as discussed above, see also 

0H. 

In the left panel of Fig. we compare the two approximate expressions (|75fl and l|8U|l for the averaged torque with 
the results obtained from the numerical solution of the FPE. Displayed is the dependence of the torque on the noise 
intensity D which is different in the two approximations. In order to keep the expansion parameter e ~ H/D roughly 
constant when changing D we have scaled the magnetic fields strengths with D. As can be seen, up to values of 
e ~ 0.5 the accuracy of l|8U|l is rather good whereas l|75|) describes the torque only qualitatively. 

Proportional change of magnetic field strength and noise intensity is equivalent to a change of the frequency u) of 
the external field in 0J. Fig. [S] hence also demonstrates the resonance- like character of the investigated ratchet effect. 
(N z ) is largest when the deterministic time scale of the external driving matches an intrinsic stochastic time scale of 
the system related to the Brownian relaxation time tb defined in (|13|l . 

The two approximate expressions for the torque (|75|l and I80|) differ in their dependence on the phase angle S. 
When using l|8U|l to fit the experimental results reported in |lfjj| we find for the Brownian relaxation time the value 
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tb — 6.4 10 _4 s instead of tb — 1.8 10 _3 s as obtained on the basis of (|75|1 in ^(j- The experiments were done with 
a ferrofiuid with 77 ~ 0.1 Pas. Using (|13[) the two results for tb hence translate in fits for the particle diameter of 
d ~ 26nm and d ~ 36nm, respectively. These values exceed the typical diameter of roughly 10 nm by factors between 
2 and 3, respectively. The main reason for this discrepancy is probably the polydispersedness of real ferrofluids, 
having a particle size distribution with a long tail (see e.g. |3j, ch.2). Hence, a whole spectrum of relaxation times 
is necessary to accurately describe the dynamics of the magnetization. However, this issue is not at the focus of the 
present investigation. We also note that the transferred torque typically increases with the particle size and hence 
the described ratchet effect is likely to be dominated by the larger grains in the ferrofiuid. 

A similar approximate calculation of the average torque is also possible for the time dependence Q ■ The result of 
the effective field approximation is 

= 45 D(Id 2 + 1) 2 ' (81) 
whereas we find from the perturbative solution of the FPE 

r(°) 



H X H^(H^>) 2 UD 2 



30 D {AD 2 + 1) 2 (36D 2 + 1) ' ^ 

The two expression only differ significantly from each other if the noise intensity D is very small. In the right panel 
of Fig. they are compared with the result from the numerical solution of the FPE l|47|) . For e ~ 0.5 the accuracy is 
again seen to be rather satisfactory 



V. CONCLUSIONS 



In the present paper we have theoretically investigated the rotational Brownian motion of colloidal ferromagnetic 
particles in an oscillating magnetic field. The central tool was the Fokker-Planck equation for the probability density 
P(e,t) of the particle orientation e at time t. Solving this equation either numerically or approximately by using 
the effective field approximation as well as a perturbative expansion we have determined the time averaged torque 
(N z ) exerted by the magnetic field on the particles. The main and a priori quite unexpected qualitative finding is 
the fact that a purely oscillating magnetic field without net rotating component can transfer angular momentum to a 
ferromagnetic grain. As basic mechanism behind this transfer, a ratchet effect was identified by which the magnetic 
field rectifies the thermal fluctuations of the particle orientation that arise due to random collisions with the molecules 
of the carrier liquid. The detailed operation of this ratchet effect in the present system was discussed on the basis of 
a weak- noise analysis of the Fokker-Planck equation for a closely related one-dimensional system in section II V Al 

Via the viscous coupling to the carrier liquid the torque on the particle is transmitted to the liquid. The combined 
microscopic torques of the huge number of individual nanoparticles then yields a macroscopic rotation of the ferrofiuid 
as a whole, as observed experimentally in 10] . In the absence of thermal noise, no average torque can arise, i. e. thermal 
fluctuations are an indispensable requirement for the rotation of the individual grains and hence of the ferrofiuid as 
a whole to occur. 

The results from the Fokker-Planck equation are completely consistent with rather general symmetry considerations 
detailed in section ITTT1 Moreover, they agree very well with simulations of the corresponding Langevin equations as 
given in |30j and quantitatively describe the experimental findings reported in |10| . 

We remark that our present system puts forward a new type of thermal ratchet device which does not fit into any 
previously known specific class of ratchet systems. Adopting the classification scheme from [l^, our present system has 
some similarity with so-called asymmetrically tilting ratchets as well as with a so-called traveling potential ratchets, 
however in the generalized sense involving two counterpropagating traveling potentials. Yet there remain significant 
differences with both these classes. In particular, we note that we are not dealing here with a periodic, asymmetric 
so-called "ratchet" -potential. Rather, at any fixed instance of time the relevant potential is perfectly symmetric about 
the instantaneous direction of the magnetic field and hence no preferential direction of rotation seems to exists. It 
is only via the time evolution that a symmetry breaking arise, which is sometimes called a dynamical symmetry 
breaking. 

In summary we hope that the present investigation has demonstrated that ferrofluids are very suitable systems 
to study various aspects of thermal ratchet behaviour, and that it may stimulate further theoretical, numerical and 
experimental work in this direction. 
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